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Abstract 

We review some approaches to the understanding of fluctuations in some 
models used to describe socio and economic systems. Our approach builds on 
the development of a simple Langevin equation that characterises stochastic 
processes. This provides a unifying approach that allows first a straight- 
forward description of the early approaches of Bachelier. We generalise the 
approach to stochastic equations that model interacting agents. Using a sim- 
ple change of variable, we show that the peer pressure model of Marsilli and 
the wealth dynamics model of Solomon are closely related. The methods 
are further shown to be consistent with a global free energy functional that 
invokes an entropy term based on the Boltzmann formula. 

A more recent approach by Michael and Johnson maximised a Tsallis 
entropy function subject to simple constraints. We show how this approach 
can be developed from an agent model where the simple Langevin process 
is now conditioned by local rather than global noise. The approach yields a 
BBGKY type hierarchy of equations for the system correlation functions. Of 
especial interest is that the results can be obtained from a new free energy 
functional similar to that mentioned above except that a Tsallis like entropy 
term replaces the Boltzmann entropy term. A mean field approximation 
yields the results of Michael and Johnson. 

We show how personal income data for Brazil, the US, Germany and the 
UK, analysed recently by Borgas can be qualitatively understood by this 
approach. 

Agent models 

We consider a simple type of agent model consisting of a set {i) of agents, 
each of which can have an intrinsic attribute, Xi. This could be a measure 
of attractiveness or repulsiveness. It could be proportional to the financial 
health or wealth of the agent. Now assume that this gives rise to some kind 
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of peer pressure or a force that tends to make 'people keep up with the 
Jones' as they say in the UK. Such an interaction force f{xi-Xj) between the 
different agents is very similar to that between molecules interact ing via an 
inter-particle potential. The strength of this kind of 'social' force, as has 
been pointed out by for example, Krawiecki, Holyst and Helbing ^ , is not 
dependent on spatial neighbourhood; rather its origin lies in communication 
between agents. The spatial topology of agents may be considered to be 
unimportant. The simplest type of peer pressure is introduced by assuming 

/ {X, - Xj) = -/ {Xj - X,) (1) 

Contrast this with the usual assumption for molecular systems where 
f{x) = f{\x\). The peer pressure force has not only a magnitude but also a 
direction depending on whether the argument, x, is positive or negative. 

Now assume the agent dynamics can be captured by the following set of 
Langevin equations: 

^=^i+j:f(x^-x,) (2) 

The stochastic terms, {^i} are defined as: 

(e. (t) it')) = 2DSijS (t - t') ; (e. (t)) = (3) 

We now compute the many agent probability distribution function 

Pr,{x,t) = (l[6[X,{t)-x,]) (4) 
A standard calculation gives: 



\i=i 



Reducing the N coordinate representation yields a hierarchy of equations 
of agent probability distribution functions: 



Krawiecki, J A Holyst and D Helbing, arXiv: cond-mat/0210044 vl 2 Oct 2002 
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9Pi {x\t) D d'^Pi {x\t) d 



-g^J dx'f{x-x')P2{x,x'\t) (6) 



dt 2 dx^ 

This hierarchy will be recognised as a dynamical generalisation of the 
BBGKY hierarchy used in the theory of liquids. Clearly a solution is only 
possible if we have a route to truncating this infinite hierarchy of coupled 
equations. The simplest way forward is to assume a mean field approxima- 
tion: 

P2{x,x'\t)^P,{x\t)P,{x'\t) (7) 
Dropping now the suffix when N=l we obtain 

dP(x\t) Dd'^P(x\t) 5 , 

The mean field force is given by: 

{x\t) = J dx'f {x - x') P {x'\t) (9) 

Now consider the agent pair force /. We have noted above (equation (1) 
that it has a directional character. Clearly a number of functions fit this 
form. A particularly simple choice is: 

/(a;)=ae-^/^ a; > , . 

f (x) = -ae^/^i x<0 ^ ^ 

If q; > 0, we see that agents with attribute Xj < Xi are attracted to agent 
i. Equally agents with attribute Xj > xi choose to avoid agent i. This mimics 
a kind of 'keeping up with Jones' social force. 

This now leads to the expression: 



oo X 

(f) {x\t) ^aj dx'e-^^'-^^/^P {x'\t) -a J d^'e-^^-^'^/^^P {x'\t) (11) 

X -oo 

This is almost, but not quite, equivalent to the form used by Marsili ^ 
who used essentially a normalised form of the above. Our form will prove 

Cccconi, M Marsili. JR Banavar & A Maritan Diffusion, peer pressure and tailed 
distributions asXiv: cond-mat/0202212v2 2 Jul 2002 
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useful in illustrating the link to the Lotka Volterra approach developed by 
Solomon and co-workers^ . 



Generalised Lotka Volterra Equations - approach of Solomon 

Using our particular choice of force defined in equation (10), we introduce 
the change of variable 

Wi = e^'/^ (12) 

Now use the theorem of Ito ^ , which states that if X is a stochastic 
process satisfying: 

dXt = (TtdZt + Htdt (13) 

Where Zt is a Wiener or Brownian process, then providing F(x) is a 
continuous, twice differentiable function Yt=F{Xt) is also a stochastic process 
and is given by 

dYt = {atF' {Xt)) dZ^ + [p.^F' {X^) + ^F" (X,) j dt (14) 
Thus the stochastic equation for our new variable Wi defined in (12) is: 

dw- y;l+C/Cl D 

Wj<Wi Wj>Wi ^ 

Note that now the stochastic function £i (t) takes only positive values 
and the additional term D may be absorbed within the mean value of this 
stochastic function. 

We refer to (15) as GLVPP (Generalised Lotka Volterra with Peer Pres- 
sure). Note that they are exact consequences of introducing the specific force 
(10) into the initial Langevin dynamics defined by equations (2). No mean 
field approximation has been used to reach this point. 



Solomon and P Richmond Power Laws of Wealth, market order and market returns 
Physica A 299 2001 188-197 

*See for example Financial Calculus: An Introduction to Derivative Pricing, M Baxter 
& A Rennie, Cambridge University Press 1996 
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For large values of Wj, the third term on the RHS of equation (15) may 
be neglected. Equally the sum over Wj < Wi may be extended over all wj. As 
a result equation (15) reduces essentially to the generahsed Lotka-Volterra 
equations proposed by Solomon for wealth dynamics. These latter equations 
have been solved using both a mean field approach and exact numerical 
methods. Both approaches show that the distribution of wealth, w, displays 
a power law or fat tail ^ . 

The comparison between our simple peer pressure model and the GLV 
model of Solomon is shown in Figure 1. The points are computed from the 
Langevin equation not the mean field theory so is an exact solution. Clearly 
the difference for the parameters used is minimal. 

The equation does not exhibit time correlations that will lead to clustered 
volatility. Solomon and Louzoun ^ have however generalised the equation es- 
sentially by replacing the random term, e (t)with the non-linear random term 
that includes a 'dynamic volatility' and proceed to show that this modifica- 
tion leads to time auto correlation functions for the volatility that exhibit 
slow long time decay in line with observations. In some respects this is similar 
to that of the Heston model 

The Heston model 

Dragulescu and Yakovenko have studied a particular stochastic model termed 
the Heston model ^ . This is characterised by two equations: 

^^fw + D{tf'we'{t) (16) 

^ = -^[D{t)-D)+^D{tf/'e^t) (17) 

The volatility, D(t), is here assumed to fluctuate about the mean, D 
with characteristic relaxation time I/7. A so-called variance noise kis also 
introduced. The second stochastic process, £^ (t) is assumed to be correlated 
with the first i.e. < e^e'^ >— p where -li nj 1. 

Malcai, O Biham P Richmond and S Solomon, Phys Rev E 66 2002 31-102 
Louzoun & S Solomon Volatility driven market in a generalised Lotka Volterra 
formalism arXiv: cond-mat/0109493 vl 26 Sept 2001 

D Dragulescu and V M Yakovenko arXiv: cond-mat/0203046 v3 5 Nov 2002 
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Known in the finance literature as the Cox-Ingersoll-Ross process, and 
in the mathematical literature as an example of a Feller process, equations 
(16) and (17) can be solved completely in terms of modified Bessel functions. 
Dragulescu and Yakovenko show that their solutions fit well Dow Jones data 
compiled over the period 1882 -2001 and 1990 - 2001. The data scales over 
time differences for the returns of 10 to 250 days for particular values of the 
parameters /, 7, D and n. The authors suggest that it will fit other 
financial time series. It remains to study the time autocorrelation functions. 

The approach has also been studied in a rather different way by Kozuki & 
Fuchikami ^ . These authors effectively solve equation (16) for the stationary 
conditional probability p{w^ D). Having also obtained the stationary solution 
to equation (17), they use standard relations to obtain the unconditional 
probability distribution, p{w) which of course agrees with the stationary so- 
lution obtained by Dragulescu and Yakovenko. 

This is an interesting approach based on empirical data for D(t) that 
suggests it may be fitted by an equation of the form (17). It extends the work 
of Bachelier but the underlying dynamics in terms of agents, as has been also 
noted by Johnson et al ^ , remains unexplained. Even more recently, Dover 
^° has shown how using a fluctuating entropy function that is then averaged 
over an ensemble of entropy sub systems, the Tsallis entropy also emerges. 

Entropy functional approach 

It is often possible for physical systems to recover the Fokker-Planck dynam- 
ics from consideration of a free energy functional for the system. Rather 
than go directly to expressions involving averages, we can go by way of Ito's 
theorem. Return to the Langevin equation (2) and consider the density 
functionpj {x, t) = 5 {x — {t)). Applying Ito's theorem yields the stochas- 
tic differential equation: 

Kozuki and N Fuchikama arXiv: cond-mat/0210090 vl 4 Oct 2002 
^NF Johnson, P JcfFcries and PMH Hui Financial Market Complexity, Oxford Univer- 
sity Press, Oxford, 2003 

^"Y Dover preprint submitted to Physica A 
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Now exploit the fact that the Dirac delta function satisfies the property: 
^^^dx^^ ~ ^^^x^"* ^^"-^ ^^^^ equation (18) reduces to equation (6) as it 

should when we average over the random variables. Using this property and 
summing over the agent index i gives a stochastic differential equation for 
the density p{x,t) — {x — (t)): 

i 

This is not yet a closed equation for the total density function p{x,t). 
However, we now follow Dean and introduce the noise term 
This has the correlation function: 

a {x, t) C {x', t')) = 2D6 {t - t') Y: Q^Q^Pi A (^', t') (19) 

Again using the property noted above of delta functions, this expression 
simplifies to: 

{i (x, t) i (x', t')) = 2D5 {t - t') ——5{x~ x') p (x, t) (20) 
We may now redefine the noise field in terms of a global noise field: 

e'(.,t) = ^^^^^^^^^ (21) 

The white noise field, -q, satisfies: 

{r]{x,t)r]{x\t')) = 2D5{x-x')5{t-t')- {r]{x,t)) = (22) 

The Gaussian noises i and i' have the same correlation functions and are 
statistically identical. Equation (??) may then be re-expressed as a closed 
Langevin equation for the global density of the system: 



(23) 

Consider now the following coarse-grained free energy functional 

F [p] = J dx J dx'p (x) V{x- x')p {x') + D J dxp (x) logp {x) (24) 



"D S Dean J Phys A Math Gen 29 No 24 21 December 1996 L613-617 
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The potential V is defined in terms of the force in the usual way, i.e./ = 
—dV/dx. The Langevin equation (23) for the density, n, is recovered now 
using the following prescription: 



dp _ d ( d SF 




(25) 



dt dx \dx dp p^p^^^^ 



Where the global noise field i' is given by equation (21). 



'Tsallis' processes 



Tsallis and colleagues have proposed the following so called entropy function 

12 . 



It is easy to see and has been pointed many times that this expression 
reduces to the Boltzmann entropy in the limit that q tends to zero although 
the fundamental basis for the postulated form of Tsallis is contentious. Nev- 
ertheless it has been used to describe a wide range of phenomena associated 
with criticality and complex systems. In particular Michael and Johnson 
have shown how this may be used to fit S&P stock market data . They 
maximise the Tsallis entropy given by equation (26) together with the fol- 
lowing constraints: 




(26) 




(27) 




(28) 



{{x 



(29) 



This yields 



^^"^'^^^Zm^^^^^^^ (g-l)[a;-a;(t)]^} 



(30) 



CJ Tsallis Stat Phys 52 1988 479 

F Michael and MD Johnson arXiv: cond-mat/0108017 vl 1 Aug 2001 
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The Lagrange multipliers Z and a are given by 



Z{t) 



^(9 - 1) P {t) 



(31) 



1 



(32) 



2a2 it) Z^~^ it) 



B(x,y) = r(x)r(y) /r(x+y) is Euler's Beta function and the usual variance 
of the distribution is 



Clearly for application to data with finite variance, the parameter q must 
lie within the range 1 to 5/3. 

This approach seems to yield some of the important stylised features 
of stock data, power laws for example. Does this approach have a link to 
Langevin processes and can we obtain a free energy functional in the manner 
of the previous section that describes the system? 

Let's return to our liquid analogue introduced at the beginning of this 
paper. In the simple approach of Einstein, it is assumed that 'collisions' of 
our colloid particle are only with the molecules that comprise the support- 
ing fluid. However imagine that the colloid particle is now actually one of 
the molecules of the fluid. So we are considering collisions between the fluid 
molecules themselves. Each molecule still undergoes random collisions how- 
ever it now seems plausible to assume that the strength of each collision is 
proportional to the density of the fluid. The simplest approximation is to 
assume an average density. However suppose that our fluid is undergoing 
strong density fluctuations of the kind one might expect close to a critical 
point. Now we might assume that the amplitude of our random collision 
process is proportional not to the global density but to some power, z, of the 
local density, n{x,i). 

With the addition of a drift force /(xj, the Langevin equation now takes 
the following form: 




q < 5/3 
q > 5/3 



(33) 
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where as before (•34') 

p{x,t) = j:s{x-Xj{t)) 
j 

The idea of using such a construction may seem a httle strange when 
thinking about hquid state theory. It will only be even approximately valid 
if the local density fluctuation has a life time much greater than the colli- 
sion process and this aspect may require more investigation. Certainly the 
process is open to criticism. However when thinking about agent models, 
such an interaction seems a plausible construction for the manner in which 
agents interact with others. The idea of local noise is a familiar one amongst 
analysts who estimate earnings and value of stock prices. Local consensus 
frequently occurs! Bear in mind that the variable, x now is not a spatial 
parameter. Rather it is an attribute of the agent. One might then follow 
the results discussed earher and choose a different form for the amplitude. 
Perhaps agents receive information not only from the 'local' community but 
also agents holding greater wealth. In this case one might choose a term of 
the kind {x, t) dx. Here we shall stick with equation (34) that at least 

has the merit, as we shall see, of leading to a particularly interesting solution. 

Using Ito's theorem, we can now deduce the dynamic equation for the 
density pf 



(Ci (t) p^ {x, t) Pi (x, t))+D— (a;, t) (x, t)) 



at ox ox ox^ 

(35) 

Taking the average to compute the probability distribution and summing 
each side over /and using the fact that the density function is a delta function 
yields: 

dP{x,t) _ df{x)P{x,t) , ^ d'P2u+l{x,t) 
~d^- +^ d^^ ^^^^ 

Once more we have a hierarchy of coupled equations. Now we may apply 
mean field theory to the second term on the RHS of (36): 

^ {Pf''^' (37) 

Introducing q = 2i+l yields 
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dP{x,t) _ difix)Pix,t)) ^ ^ d'P^x,t) 
dt dx dx'^ 

Choosing the force /~ A+Bx it is readily seen that the stationary solution 
to this equation takes the form given by equation (30). This suggests that we 
might be able to formulate a free energy functional that involves the Tsallis 
formula. Summing equation (35) over the index i gives 

dp d[f{x)p] d[ri{x,t)p^/'] ^ ^d'p^ 
dt dx dx dx"^ 

Now follow the procedure used earlier. It is straightforward to show that 
the noise 

Cix,t)^j:-^iCiit)p''ix,t)p,{x,t)) (40) 

i 

has the same statistical property as the global noise 



d 



T] {X, t) ^) 



CM^^—--^^ (41) 



dx 

Where the noise rj is defined as before through equation (22). 
So once again we have a closed Langevin equation for the evolution of the 
density function, n: 

dp_ d[f{x)p] d[vix,t)p'^/']^^d^ ^^^^ 



dt dx dx dx"^ 

Forming the average quantities yields equation (36) as expected. More 
to the point this equation also follows from the following coarse grained free 
energy defined below together with the procedure defined by equation (25) 
(with i replaced by ce) 

F[p] = J dxV (x) p (x) + [l- j dxp'i {x, t)^ (43) 

Thus we have now a free energy functional for our system where the role 
of entropy previously played by a Boltzmann expression is replaced by the 
Tsallis formula. 
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In approaching these complex systems via the Langevin equation we can 
begin to see the nature of the approximations involved. Already we have 
noted that the local noise could be characterised by a different form includ- 
ing, for example, correlations with agents of greater wealth. In the case 
of a liquid we might consider extending the influence to include molecules 
in the wider neighbourhood. For example we might replace (,T)with say 
c?a;'e~"'^~^''p'^ (a;')where a is a parameter characterising the range of in- 
fluence. Clearly a = yields the Boltzmann hmit whereas a ^? is consistent 
with the Tsallis limit. Similarly for socio economic systems one plausible 
approximation could be to replace {x) with/^ rfx'e""*^^"^ ^ p^ (a;')which as- 
sumes a peer pressure where only 'wealthier' people have an influence on 
individual behaviour. The reader may think of other variants. The point 
we make is to suggest that the Tsallis expression is not a unique expression. 
Rather it may be an approximation albeit a rather clever one that accounts 
in a subtle manner for important correlations in these complex systems. 



Wealth distributions 

In a recent paper Borgas studied data for wealth distributions in the USA 
(2002), Brazil (1996), Germany (1998) and the UK (1998). The cumulative 
distributions, he found could be fltted well with a q q' Gaussian obeying the 
differential equation: 



1 dP 

= - W- P')P'> - P'P"'; l<q'<q and < < /3 (44) 

2x dx 

Full details are given in Borgas' paper. Here we note the values he ob- 
tained for the exponents q and q': 



(45) 





q 


q' 


USA 


3.8 


1.7 


Brazil 


3.5 


2.1 


Germany 


2.9 


1.0 


UK 


3.1 


1.0 
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E P Borges arXiv: cond-mat/0205520 v2 24 May 2002 
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Empirical data suggests that q' remains constant whereas over the past 
30 years there has been a small increase in the values for q. 

Return to equation (38). Suppose the force acting on agent i arises from 
not only pair forces but also higher order many body forces. In this case it 
can be expressed in terms of a many body expansion, 

f{x,[p]) = a + bp + cp^ + .... (46) 

And a is identified with external effects; b with one body interactions, 
c with 2 body interactions, etc. If we invoke the mechanism used in the 
Sznajd model, then we assume agents do not respond to single body forces. 
Only two body forces affect the agent dynamics. So we must set b=0. Now 
examine the stationary solution of equation (38): 

^ = aP + cP^ (47) 
ax 

If q ~ 1 then the exponents are close to those of the European countries, 
UK and Germany. If q j 1 we may write the term on the LHS of equation 
(47) as qP'^~^dP/dx and obtain 

HP 

^ ^ a'P^~i + c'P^-'i (48) 

dx 

A value of q ~ 0.3 now yields exponents of 1.7 and 3.7 which within 
experimental error seem to be compatible with the results for the US. More 
detailed analysis of the raw data is required to see how close equation (48) 
really is for all the data including Brazil. Furthermore, the so-called 'an- 
kle' that is evident in the very high-income data of Borgas remains to be 
explained. However this approach suggests that income data is governed by 
similar mechanisms to investments and trading. Furthermore agents within 
the American continent, are much more responsive to local noise and rumour 
than agents in Europe. Job switching to secure higher incomes in the US is 
known to be more dynamic than in Europe. Certainly that great American 
investor the so-called 'Sage of Omaha', Warren Buffet is on record as saying 
he lives in Omaha simply to avoid being influenced by rumours in Wall St. 
Perhaps Europe has traditionally been less influenced by local rumour and 
noise. More detailed studies might cast more light on this feature. 
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Figure caption 
Figure 1. 

Normalised cumulative distribution functions for the wealth in both the 
standard GLV model of Solomon and the LV model with peer pressure (equa- 
tion (15)). The fit for the highlighted region (in yellow) is in good agreement 
with the analytical (mean-field) results. The last part of the distribution dis- 
plays a sharper drop-off. Finite-size effects may account for the discrepancy 
between analytical and numerical results. 
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